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Abstract. New experiments for water at the surface of proteins at very low 
temperature display intriguing dynamic behaviors. The extreme conditions of these 
experiments make it difficult to explore the wide range of thermodynamic state points 
needed to offer a suitable interpretation. Detailed simulations suffer the same problem, 
where equilibration times at low temperature become unreasonably long. We show how 
Monte Carlo simulations and mean field calculations of a tractable model of water help 
interpret the experimental results. Here we summarize the results for bulk water and 
investigate the thermodynamic and dynamic properties of supercooled water at an 
interface. 
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1. Introduction 

Water is essential in biology, because it participates in nearly every process necessary for 
life (including cell metabolism, transport of nutrients and residues, protein conformation 
changes, etc.), and is the most common solvent in chemistry. It regulates a large variety 
of processes, including atmospheric phenomena, the formation of geophysical structures, 
the propagation of cracks in stones and cement, the sliding of glaciers, the transport 
in plants, and is ubiquitous in the universe as ice in the interstellar space. In all these 
examples the properties of water are essential to understand what is observed. 

Nevertheless, water has proven to be a complex puzzle to many researchers 
for its anomalous thermodynamic and dynamic properties at room temperature and 
pressure. For example, by decreasing temperature T at pressure P = 1 atm, water's 
volume fluctuations, proportional to the isothermal compressibility Kx, increase below 
T = 46°C, and entropy fluctuations, proportional to the isobaric specific heat Cp 
increase below T = 35°C, while in normal liquid any fluctuation decreases when T is 
decreased [HO [3]. These water's anomalies grow upon cooling and increase in number. 
For example below T = 4°C the cross-fluctuation of volume and entropy, proportional 
to the isobaric thermal expansion coefficient ap, becomes negative [1], while it is always 
positive in normal liquids where the entropy decreases when the volume decreases [5]. 

By decreasing T even more, it is experimentally possible to supercool liquid water 
down to Th = — 41°C at 1 atm and to Tjj = — 92°C at 2000 atm, where the liquid is 
metastable with respect to crystal phases [lj. These extreme conditions are not unusual 
in nature, where water exists in its liquid form at — 20°C in insects, — 37°C in clouds or 
— 47°C in plants [6]. Below Th the homogeneous nucleation of crystal ice occurs in a 
time too short to allow any measurements. But even the crystal phase of water is not 
simple. In fact, water is a polymorph with at least sixteen forms of crystal ices, the last 
one, Ice XV, was discovered in 2009 [7j. 

However, at very low T, crystal water is not the only possible kind of ice. By 
rapidly quenching liquid water below — 123°C |8j, or by condensing the vapour at low T 
[9], or by compressing crystal ice at low T [10], or by irradiation (with ions for example 
it is possible to solidify water as an amorphous, or glass, i. e. a form that has 
the elastic properties of a solid, but the structure of a liquid with no long-range order 
[T2] . As for the crystal state, the amorphous state of water is also not unique. Water 
is a polyamorph with at least three different amorphous states: low-density amorphous 
(LDA), discovered in 1935 [13], high-density amorphous (HDA), discovered in 1984 [14J, 
and very high-density amorphous ice (VHDA), discovered in 2001 [15] . 

All these anomalies of water are a consequence of the properties of the hydrogen 
bond network that they form. The hydrogen bond interaction is characterized by a 
preferred geometrical configuration, that at low T and P is approximately a tetrahedron 
of four molecules around a central one, with an angle varying around 106.6° (slightly 
smaller than a tetrahedral angle of 109.47°) and a distance oscillating around 2.82 A[5]. 
The local arrangement, including the number of nearest neighbours, can change with T 
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and P. In particular, in 2000 Soper and Ricci observed at 268 K, compressing from 26 to 
400 MPa, a continuous transformation from low-density liquid (LDL) local arrangement 
of water with an open, hydrogen-bonded tetrahedral structure, to high-density liquid 
(HDL) local arrangement with nontetrahedral 0-0-0 angles and a collapsed second 
coordination shell with broken hydrogen bonds, and a change in density of about 73% 
PJ- 

1.1. Thermodynamic Interpretations of Water Behavior 

All the above results are consistent with theories that propose different mechanisms and 
different phase behaviors in the supercooled region. They can be summarized in four 
possible scenarios for the P — T phase diagram. 

(i) In the stability limit (SL) scenario [17] the behavior of the superheated liquid 
spinodal, i. e. the limit of stability of the liquid with respect to the gas, and the stretched 
water, i. e. water under tension as in a plant fibers, are related. In particular, it is 
hypothesized that the limits of stability of these two regions are continuously connected 
at negative pressure, forming a re-entrant curve toward the positive P region below 
T H (P). The response functions, including K T , Cp and ap, diverge when T is decreased 
a positive P as a consequence of the approaching of the re-entrant spinodal line. 

(ii) In the liquid-liquid critical point (LLCP) scenario [18] it is hypothesized the 
existence of a LDL-HDL first-order phase transition line with negative slope in the 
P — T plane and terminating in a critical point C Below the critical pressure Pc 
the response functions increase on approaching the Widom line (the locus of correlation 
length maxima emanating from C into the one-phase region) [5], and for P > Pc by 
approaching the HDL-to-LDL spinodal line. The possibility with Pc < have also 
been proposed [19J. 

(iii) In the singularity-free (SF) scenario [20J it is hypothesized that the low-T 
anticorrelation between volume and entropy gives rise to response functions that increase 
upon cooling and display maxima at non-zero T, but do not display any singular 
behavior. Specifically, Sastry et al. [20] show that this is a direct consequence of the 
fact that water's line of temperatures of maximum density (TMD) has a negative slope 
in the (T, P) plane. 

(iv) In the critical-point free (CPF) scenario [21J it is hypothesized that a first- 
order phase transition line separates two liquid phases and extends to P < toward 
the superheated limit of stability of liquid water. This scenario effectively predicts a 
continuous locus of stability limit spanning the superheated, stretched and supercooled 
state, because the spinodal associated with the first-order transition will intersect the 
liquid-gas spinodal at negative pressure. No critical point is present in this scenario. 

Since all these scenarios are consistent with the available experimental data, a 
natural question is if we can design an experiment that would discriminate among 
them. Unfortunately, many scientists have discovered that an answer to this question is 
a difficult challenge [21] . In fact, experiments on bulk water are hampered by freezing 
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below Th, and no measurements on bulk liquid water can be performed with our present 
technology below this temperature. Nevertheless, the different proposed theories have 
different implications on phenomena such as the cold denaturation (and stop of activity) 
of proteins at low T, an important issue in cryopreservation, cryonics, cryostasis and 
cryobiology. 

1.2. Interfacial Water 

One possible strategy to probe supercooled water at very low T is to consider 
water at an interface. Water adsorbed onto the surface of proteins or confined in 
nanopores freezes at much lower T than bulk water, giving access to a low-temperature 
region where interfacial water is still liquid, while bulk water would not be [22J. 
In many cases of interest for practical purposes in biology, geology or industrial 
applications, water is hydrating a surface or is confined. As a consequence, fundamental 
research in physics and chemistry has been performed in recent years with experiments 
[HI [231 [2H [251 1261 [271 12HI [211 [221 ED], theories and simulations [3D E21 EH E3] . 

During the last years experiments on water in Vycor micropores [35], in nanopores 
of MCM-41 silica [36J, EH [38] , of sol-gel silica glass [39], of NaA zeolites [3D], or of double- 
wall carbon nanotubes [41] have contributed to the investigation of water dynamics in 
confinement. In particular, confinement in hydrophilic MCM-41 silica nanopores of 1.8 
and 1.4 nm diameter allows to study water dynamics down to 200 K where quasielastic 
neutron scattering reveals a crossover at T w 225 K in the average translational 
relaxation time from a non-Arrhenius behavior at high T to an Arrhenius behavior 
at low T [32] • A similar crossover is also observed for the self- diffusion coefficient of 
water by nuclear magnetic resonance at T ^ 223 K [43J . By increasing from 400 bars to 
1600 bars the external pressure applied on a sample of MCM-41 silica nanopores with 
1.4 nm diameter at full hydration level of 0.5 g of H2O per g of silica, it has been observed 
that the crossover occurs at lower T, reaching T « 200 K at P = 1600 ± 400 bars and 
disappears at higher P |44j. 

Quasielastic neutron scattering experiments show the same crossover for the average 
translational relaxation time of at least three different systems: (i) water hydrating 
lysozymes, at hydration level h = 0.3 g of H2O per g of dry lysozyme, for T w 220 K, 
a temperature below which the protein has a glassy dynamics [30], (ii) DNA hydration 
water, at hydration level corresponding to about 15 water molecules per base pairs, for 
T m 222 K, at which DNA displays the onset of anharmonic molecular motion |45j, 
and (iii) RNA hydration water, at a similar hydration level, for T ph 220 K, where 
both RNA and its hydration water exhibit a sharp change in slope for the mean-square 
displacements of the hydrogen atoms [46J, 

All these results can be interpreted as a consequence of a structural rearrangement 
of water molecules associated with a LDL-HDL critical point [33]- In fact, along the 
Widom line in the supercritical region of the LDL-HDL critical point, the changes in 
the hydrogen bond network are consistent with the dynamic behavior observed in the 
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experiments, as we will discuss in the following sections. 

This interpretation has been criticized [47] on the base of similar crossover observed 
for water confined in molecular sieves [48] or for water mixtures [49] and water solutions 
[50] . It has also been proposed a possible interpretation as a consequence of the dynamics 
of (Bjerrum-type) defects due to orientationally disordered water molecules that are 
hydrogen bonded to less than four other water molecules [51], [52] . 

2 H-NMR studies on hydrated proteins, at a comparable hydration level as in [44j . 
show no evidence for the crossover at 220 K and indicate that water performs thermally 
activated and distorted tetrahedral jumps at T < 200 K, which may be related to a 
universal defect diffusion [2H]. Also, dielectric spectroscopy studies of hydrated protein 
show a smooth temperature variations of conductivity at 220 K and ascribe the crossover 
observed in neutron scattering to a secondary relaxation that splits from the main 
structural relaxation [291 E3] . 

On the other hand, numerical simulations for bulk water show that crossing the 
Widom line emanating from a LDL-HDL critical point, the structural change in water is 
maximum, as emphasized by the maximum in specific heat, and the diffusion constant 
has a crossover [54] • This result is observed also in water hydrating lysozyme or DNA, 
where the dynamic transition of the macromolecules occurs at the temperature of 
dynamic crossover in the diffusivity of hydration water and also coincides with the 
maximum change in water structure [55] . 

A crossover from high-T non-Arrhenius to low-T Arrhenius behavior is observed 
also in simulation of water hydrating lysozyme powder in the translational correlation 
time and in the inverse of the self-diffusion constant, in agreement with the neutron 
scattering experiments, at about 223 K [56]. The activation energy for the Arrhenius 
regime is found to be of about 0.15 eV, as in the neutron scattering experiments [56J. 
Also, simulations of water hydrating elastin-like and collagen-like peptides show this 
crossover, but with a weaker change in the slope and an Arrhenius activation energy of 
about 0.43 eV, consistent with dielectric spectroscopy and nuclear magnetic resonance 
studies (57J EE]. 

It is therefore difficult not only in the experiments, but also in the models to get 
a clear answer about the relevant dynamic mechanisms and their relation with the 
thermodynamics in water at interfaces. Moreover, the relation between confined water 
and bulk water remains not fully understood. For this reason models that are tractable 
with a theoretical approach are particularly appealing in this context. With these 
models, in fact, simulations can be compared with analytic calculations to develop a 
consistent theory. 

2. Cooperative Cell Model for a Monolayer of Water 

We consider the case of water in two dimensions (2D). This case can be considered as an 
extreme confinement of one single layer of water between two repulsive (hydrophobic) 
walls when the distance between the walls is such to inhibit the crystal formation [58J. 
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In fact, it has been shown that the relevant parameter to avoid the transition to a 
crystal phase is the distance between the confining wall and not the characteristics of 
the hydrophobic interaction with the wall [58J. 

Another case in which the study of a monolayer of water is relevant is when a 
substrate of protein powder is, on average, hydrated only by a single layer of water, 
and the proteins do not undergo any configurational transformation and/or large scale 
motion [59J. In these conditions, for a hydrophilic protein surface, we can assume that 
the effect of the water-protein interaction is to attract water on a surface that, by 
constraining the water molecule positions, inhibits its crystallization. 

A very desirable feature of a model for a liquid is transferability. The parameters 
and effective interactions of a model are optimized to precisely reproduce static and 
dynamic properties of the liquid at one particular thermodynamic state point. The 
quality of the model is measured by the range of validity of its predictions in other 
state points. Unfortunately, there is no water model that is truly transferable, nor can 
reproduce all the properties of water [60]. Many routes have been explored to solve this 
issue. Molecular polarizability [6TJ [62] is one way to introduce effects not considered by 
standard pairwise additive potentials. However, polarizable models are computationally 
very expensive and provide only a partial solution [63] . An alternative way is to include 
many-body effects into the potential. In the following we define a model with an effective 
many-body interaction introduced through a cooperative hydrogen bond term. 

2.1. Definition of the Model with Cooperative Interaction 

We consider N molecules in a volume V with periodic boundary conditions (p.b.c.) in 
two dimensions, and the size of about one single molecule (and no p.b.c.) in the third 
dimension. We initially consider the case in which the molecules are distributed in a 
homogeneous way, with each molecule i G [1, N] occupying the same volume V/N larger 
than a hard-core volume Vq ~ 102 A 3 due to short-range electron clouds repulsion. 
We consider the case in which each molecule has coordination number four, consistent 
with the tendency of a water molecule to minimize its energy by forming four hydrogen 
bonds. 

The interaction Hamiltonian for water molecules is [6U [65], EU [5] 

<i,i> (h,l)i 

where Uo(r) denotes the sum of all the isotropic interactions (e. g. van der Waals) 
between molecules at distance r = (V/N) 1 / d , represented by a Lennard- Jones potential 
with attractive energy e plus a hard-core at distance r = (fo) 1 ■ 

The second term (with 5 a ,b = 1 if a = 6 and 5 at b = otherwise, and denoting 
that i and j are nearest-neighbors) accounts for the directional contribution to the 
hydrogen bond energy with strength J, where cr^- = 1, . . . , q is a (Potts) variable 
representing the orientational state of the hydrogen (or the lone e~) of molecule i facing 
the lone e~ (or the hydrogen, respectively) of the molecules j. For the sake of simplicity 
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we do not distinguish between hydrogen and lone e~, associating to each molecule four 
equivalent bond indices Cy. We choose the parameter q by selecting 30° as the maximum 
deviation from a linear bond, i. e. the O — H....0 angle is less than 30°, as estimated 
from Debye- Waller factors [67J EE]. Hence, q = 180°/30° = 6 and every molecule has 
q A = 6 4 = 1296 possible orientations. The effect of choosing a different value for q has 
been analyzed in [5]. 

The third term (with (k, l)i indicating each of the six different pairs of the four 
bond indices of molecule i) represents an interaction accounting for the hydrogen bond 
cooperativity and giving rise to the 0-0-0 correlation [25] , locally driving the molecules 
toward an ordered (tetrahedral in the bulk) configuration with lower energy. 

By defining the energy per hydrogen bond (between cy and o~ji) as the sum of 
the interactions in which two bonded molecules (i and j) are participating, we obtain 
Ebb = e + J + mJ a /2, where m = 0, . . . , 6 is the number of cooperative interactions 
in which that bond variables (cry and o~ji) are partaking. If we choose as parameters 
e = 5.8 kJ/mol (consistent with the value 5.5 kJ/mol of the estimate of the van der 
Waals attraction based on isoelectronic molecules at optimal separation [69]), J = 2.9 
kJ/mol and J a = 0.29 kJ/mol, the values of Ebb ranges between 8.70 and 9.6 kJ/mol 
depending on m. However, a definition of E'hb based on a cluster of 5 or 8 bonded 
molecules in d =3-dimensions increases this range up to 17 or 18 kJ/mol, respectively. 
Therefore, Ebb depends on the environment (the value of m and the number of molecules 
bonded in a cluster), as observed in computer simulation of the crystalline phases of ice 
[70], and has values within the range 6.3 [71] — 23.3 kJ/mol [75] . proposed on the base 
of experiments. 

Experiments show that formation of the hydrogen bonds leads to an open — locally 
tetrahedral — structure that induces an increase of volume per molecule [731 EE] ■ This 
effect is incorporated in the model by considering that a full bonded molecule, i. e. a 
molecule with four hydrogen bonds, has a molecular volume larger than a non-bonded 
molecule by an amount 

Av = 4vbb, (2) 

where vbb is the volume increase per H bond. Hence, if 

Nbb = riiUjS^^ (3) 
<i,i> 

is the total number of hydrogen bonds in the system, the hydrogen bond contribution 
to the total volume is 

AV = iW HB . (4) 

We adopt r = 2.9 A consistent with the expected value of the van der Waals radius [74J, 
and vbb = 0.5t> , with v = Tq, corresponding to a maximum hydrogen bond distance 
of about 3.3 A, consistent with the range of a water molecule's first coordination shell, 
3.5 A, as determined from the oxygen-oxygen radial distribution function [75j . 
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3. The Phase Diagram 

The model is studied using both mean-field (MF) analysis and Monte Carlo (MC) 
simulations. The MF approach has been describe in details in Ref.s [5], [76]. It consists 
of expressing the molar Gibbs free energy in terms of an exact partition function for a 
portion of the system made of a treatable number of degrees of freedom. We take into 
account the effect of all the rest of the system as a mean field acting on the border of 
this portion [HiESEgEElEIJIi^^ 

MC simulations are performed at constant N, P, T, allowing the volume Vmc of 
the system to fluctuate as a stochastic variable. The average distance r between the 
molecules is then calculated as r/r = (Vmc/vqN) 1 ^, where d = 2 in two dimensions. 
The total volume of the system is by definition 

V = V MC + AV, (5) 

where AV is the hydrogen bond volume contribution in Eq.Q. Note that AV is not 
included in the calculation of r to avoid MF-type long-range correlation in volume 
fluctuations in the MC simulations. 

For the parameters choice J/e — 0.5, J a /e = 0.05 and uhbA'o — 0.5, we find that 
the density p as a function of T at constant P displays a discontinuous change at high 
T and low P corresponding to the gas-liquid first order phase transition anding in a 
critical point where the discontinuity disappears (Fig. [Th). By decreasing T, the density 
reaches a maximum, that in real water at atmospheric pressure occurs at 4°C. At lower 
T, in the supercooled state, and higher P we find another discontinuity in density, this 
time with a lower density at lower T (Fig. [lb). The system at these supercooled T 
displays a first order phase transition from HDL to LDL, as hypothesized in the LLCP 
scenario (Fig. 2). 

3.1. Effect of Hydrogen Bond Cooperativity on the Behavior of Water 

The experiments for confined water have boosted the debate over the supercooled phase 
diagram of water, motivating the proposal of the CPF scenario hypothesized by C. A. 
Angell [21], as described above. This new scenario leads to questions such as 

(i) How to understand the new Angell hypothesis? 

(ii) How to connect it to the other three existing hypotheses? 

A recent work by Stokely et al. [83] succeeds in answering both questions (i) and 
(ii). Specifically, it is shown that all four existing hypotheses are cases of the cooperative 
water model. Thus no matter which hypothesis may be correct (if any is correct), it 
is possible that the underlying mechanism is basically the same — the thermodynamic 
properties of water can be accounted for by considering two main contributions to the 
hydrogen bond interaction: (a) the directional (or covalent) contribution (parametrized 
by J in the model) and (b) the three-body (or cooperative) contribution (parametrized 
by J rj in the model). 



Phase Transitions and Dynamics in Bulk and Interfacial Water 



9 




Figure 1. The density p (in units of l/«o) as a function of the temperature T (in 
units of e/ks where ks is the Boltzmann constant) for different values of pressure P 
(in units of e/vo) as calculated from MC simulation of a system with N — 15625 water 
molecules. The parameters of the model are J/e = 0.5, J a jt — 0.05 and wr-b/vo ~ 0.5. 
(a) At high T and for (from bottom to top) values of Pvo/e from 0.1 to 0.9, we observe 
for Pno/e < 0.8 a discontinuity in the density corresponding to the first-order gas- 
liquid phase transition, with a critical P at about (0.75 ± 0.05)e/i>o and critical T at 
about (2.2±0.1)e/&B. Note that if we choose as model parameters e = 2.5 kJ/mol and 
r = 3.2 A, we get an estimate P c > = 22.7 ± 1.5 MPa and T c , = 661 ± 30 K consistent 
with the real water critical point at about 22.064 MPa and 647 K. (b) At low T and 
for (from bottom to top) Pvq/c from 0.5 to 0.9, for Pvo/e > 0.55 a discontinuity in p 
marks the first-order LDL-HDL phase transition. 
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Figure 2. The finite size behavior of the maximum of compressibility K^ ax as a 
function of the number of water molecules N for pressure P = 0.8e/«o at low T shows 
a linear increase as expected at a first order phase transition. 
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By MF calculations and MC simulations, Stokely et al. [S3] demonstrate that the 
balance between contributions (a) and (b) determines which of the four hypotheses 
presented in section 1.1. best describes experimental facts. Since the characteristic 
energy associated with these two contributions can be estimated, the work allows to 
begin to validate or contradict each hypothesis on an experimental basis. 

Specifically, by fixing the parameters J/e — 0.5 and vhb/^o = 0.5, and varying 
the parameter J a /e, it is possible to observe that the cooperative model reproduces 
all four scenarios of section 1.1. The overall picture that emerges is one in which the 
amount of cooperativity among H bonds (J a /e), in relation to the H bond directional 
strength (J/e), governs the location of a LLCP, hence which scenario is realized. 
For zero cooperativity, the temperature Tc where K™ x and ctp ax diverge is at zero 
temperature, and no liquid-liquid transition exists for T > 0-the SF scenario. For very 
large cooperativity, C lies outside the region of stable liquid states, and a liquid-liquid 
transition extends to the entire (supercooled and superheated) liquid phase-the CPF/SL 
scenario. For intermediate values of H bond cooperativity, Tc varies in a smooth way 
between these two extremes-the LLCP scenario. Due to the anticorrelation between 
the volume and entropy associated to the H bonds, the larger Tc-, the smaller Pc, 
eventually with Pc < for larger cooperativity. These cases are summarized in Fig. [3j 

4. Water at Interfaces 

To elucidate the relation between the protein dynamic crossover at about 220 K and 
the dynamic crossover observed for the average translational correlation time in the first 
layer of protein hydration water [301 H5j HQ HI] , we perform MF calculations and MC 
simulations of the cooperative model of water of section 2.1. Since we are interested in 
cases at low humidity, we consider the case of a water monolayer hydrating an immobile 
surface of globular protein that, forcing the water molecules out of place with respect 
to crystal configurations, inhibits the crystallization. We focus on the hydrogen bonds 
dynamics, regardless if the hydrogen bonds are formed with the protein or among water 
molecules. 

4-1. The Hydrogen Bond Dynamics for Hydrated Proteins 

Following the work of Kumar et al. [80], [791 EE] , Mazza et al. [59J study the orientational 
correlation time r associated with the hydrogen bond dynamics of the model is in section 
2.1. They confirm the occurrence of a dynamic crossover at a temperature of about 
ksT '/e ~ 0.32 corresponding to the T of maximum variation of the number of hydrogen 
bonds, that in turn corresponds to the Widom line. They also confirm that the crossover 
is from a non-Arrhenius behavior at high T to another non-Arrhenius behavior [82J, 
that closely resembles an Arrhenius behavior around the crossover. 

These results are consistent with those from simulations of other models for 
hydrated proteins where a crossover in the translational dynamics is observed [55] . The 
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Figure 3. The pressure P vs temperature T phase diagram of the cooperative 
water model for different values of J a /e. (a) Mean field results showing the low T 
phase diagram with the LDL-HDL phase transition lines (solid lines, where Kt is 
discontinuous) and Widom lines (dashed lines, where Kt has a finite maximum) for 
varying J a /e from (rightmost) 0.04, 0.03, 0.02, 0.01, 0.005 and (leftmost). For each 
value of J a /e solid circles indicate the LDL-HDL critical point C where the response 
function (such as Kt) diverge. Hashed circle indicates the state point at T — where 
Kt diverges when J a /e = 0. (b) The MC phase diagram for varying J a je for N = 10 4 
water molecules. At high T the system displays a liquid-gas first-order phase transition 
line (continuous line with full circles) ending in a liquid-gas critical point C (full circle), 
from which departs the liquid-to-gas spinodal line (dashed line with open diamonds). 
At lower T, the retracing line with open squares marks the temperatures of maximum 
density (TMD) along isobars. All these loci do not change in an appreciable way with 
the value of J a /e. The phase diagram at lower T, instead, show a strong dependence 
from Jc/e. For J a /e = 0.5, we find for any P above the spinodal line a first-order phase 
transition line between a HDL (at high P) and a LDL (at low P) phase (continuous 
line with open circles). This is the CPF scenario [21]. The analysis of the HDL-to-LDL 
(not shown) reveals that the limit of stability of the HDL phase retraces to positive 
P as in the SF scenario [17] • For J a /e = 0.3 we observe that the HDL-LDL phase 
transition ends to a critical point C (continuous line with open circles ending in a full 
large circle) at negative P, as in the LLCP scenario suggested in Ref. [19J. From C 
a Widom line (dashed line with full circles) departs. For J„/e — 0.05, C occurs at 
positive P, has hypothesized in Ref. [H] and the Widom line (dashed line with full 
squares) extends to lower P. For J a /e = 0.02, C approaches the T = axis, as well 
as the Widom line (dashed line with full diamonds), going toward the limit of the SF 
scenario [213] , 

difference here is that (i) the crossover is for the dynamics of the hydrogen bonds, (ii) 
in the cooperative water model the crossover can be calculated from MF and an exact 
relation can be found between the crossover and the Widom line, and (iii) the model can 
test different hypotheses. In particular, Kumar et al. have shown that the crossover at 
ksT/e ~ 0.32 is independent of whether water at very low temperature is characterized 
by a LLCP or is SF. In fact, the crossover is a consequence of the sharp change in the 
average number of hydrogen bonds at the temperature of the specific heat maximum, 
that occurs in both scenarios. Kumar et al. were able also to make predictions about 
the P-dependence of quantities characterizing the crossover at A^T/e ~ 0.32: (i) the 
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time scale of the crossover, showing that it is independent of P (isochronic crossover); 

(ii) the activation energy of the apparent Arrhenius behavior at low T and (iii) the 
crossover temperature, showing that both (ii) and (iii) decrease linearly upon increasing 
P [80]. These predictions have been confirmed by Chu et al. [EH [82] in a study on the 
dynamics of a hydrated protein under moderately high pressures at low temperatures 
using the quasielastic neutron scattering method. They relate these predictions (i)- 

(iii) to the mechanical response of the protein to an external force, that is the average 
elastic constant calculated from the mean square displacement of the protein atoms. In 
particular, the degree of "softness" of the protein, related to the enzymatic activity, is 
preserved at lower T if the pressure is increased [HI]. However, a criterion proposed 
in Ref. [80J for discriminating which scenario better describe water on the basis of the 
crossover at ksT/e ~ 0.32, cannot be tested in the experiments since the predicted 
difference between the two scenarios (of the order of 1%) is within the error bars of the 
measurements [82J. 

The answer to the puzzle of which of the scenarios better describe water might be 
related to the very recent experimental discovery of another crossover for the hydrogen 
bond r to an Arrhenius behavior at very low T, of the order of 180 K at hydration 
h = 0.3 g H 2 0/g. This crossover has been observed by Mazza et al. at ksT/e ~ 0.07 
[59], in relation to an ordering process of the hydrogen bonds leading to the HDL-LDL 
critical point. The study has been possible thanks to the use of a highly efficient cluster 
MC dynamics [85l [86] . This very-low T crossover would reduce even more the T at 
which the proteins preserve their "softness" , essential for their correct functionality. 

4-2. Water Monolayer in Hydrophobic Confinement 

By considering partially hydrated hydrophobic plates at a distance such to inhibit the 
crystallization of water at low T, Franzese and de los Santos [HZ] have show that water 
has a glassy behavior [HE] for both the translational and rotational degrees of freedom 
when cooled down to a low P. This result is consistent with simulations of TIP4P 
water forming a quasi-2d amorphous when confined in a hydrophobic slit pore with 
wall-to-wall separation just enough to accommodate two molecular layers [89] . 

At higher P the hydrogen bond network builds up in a less gradual way, allowing the 
system to equilibrate the rotational degrees of freedom also at very low T, but not the 
translational degrees of freedom. This effect is emphasized by the appearance of many 
dehydrated regions [87], as also observed in water confined between two protein-like 
hydrophobic flattened surfaces at distances ranging from 0.4 to 1.6 nm [33J. 

When P is close to the LLCP value, the cooperativity of the hydrogen bond network 
induces a strong non-exponential behavior [90J for the hydrogen bond correlation 
function. However, both rotational and translational degrees of freedom equilibrate 
within the simulation time. At higher P the rotational correlation function recovers the 
exponential behavior and the diffusion of the system allows the formation a large dry 
cavity, while the rest of the surface is well hydrated. It is interesting to observe that the 
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cooperative model allows to calculate, in the MF approximation, the diffusion constant 
at any T and P [91]. 

The hydrophobic confinement has effects also on water thermodynamics. It shifts 
the HDL-LDL phase transition to lower temperature and lower pressure compared to 
bulk water when the confinement is between plates [92, 93J. Moreover, it shifts both 
the line of maximum density and the liquid-to-gas spinodal toward higher pressures 
and lower temperatures with respect to bulk when the confinement is in a hydrophobic 
disordered matrix of soft spheres [33] • This result is confirmed also in the analysis 
performed by using the cooperative water model in confinement between hydrophobic 
hard spheres [93] • However, the effect of the matrix on the HDL-LDL critical point is 
less clear and is presently under investigation. 

5. Conclusions 

The effect of confinement is of great interest to biology, chemistry, and engineering, 
yet the recent experimental and simulations results are object of an intense debate. A 
better understanding of the physico-chemical properties of liquid water at interfaces 
is important to provide accurate predictions of the behavior of biological molecules 
[95] . including the folding-unfolding transitions seen in proteins [%l [97J [98], and the 
dynamical behavior of DNA [33]. However, it is still unclear whether such behaviors 
are inherent in the structure of such molecules, or an effect of water in which they are 
always found, or due to the interactions between the two. 

To get insight into this subject the formulation of a model that allows the 
development of a theory could be useful to find functional relations connecting different 
observables. The advantage of this approach is to have two independent ways of 
approaching the problem, one theoretical and the other numerical. 

We have presented here several recent results obtained with a cooperative water 
model suitable for studies with mean field theories and with N, P, T simulations with 
thousands of molecules. The model has been studies in the context of water monolayers 
on hydrated proteins, between hydrophobic surfaces or in a hydrophobic matrix. 

Some of the conclusions reached with this model are the following. 

• The different scenarios proposed to interpret the low-T behavior of water are 
instances of the same mechanism, with different values of the directional (covalent) 
strength and the cooperative (many-body) interaction of the hydrogen bonds. The 
parameters that can be estimated from the experiments suggest that the scenario 
with the LLCP is the most plausible for water. 

• Previous experiments showing one dynamic crossover in the water monolayer 
hydrating proteins, RNA and DNA are consistent with (at least) two scenarios. 

• The possibility of a second dynamic crossovers detectable at lower T and lower 
hydration level would be consistent only with the LLCP scenario, because its origin 
would be related to the ordering of the hydrogen bond network. 
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• A consequence of the occurrence of a LLCP should be detectable when the 
translational and rotational dynamics of water are studied for a monolayer in a 
hydrophobic confinement. In particular, the rotational dynamics should appear 
with the strongest non-exponential behavior in the vicinity of the LLCP, as an 
effect of the cooperativity. Moreover, the slow increase of the number of hydrogen 
bonds at ow T and low P is the cause of the formation of an amorphous glassy state 
when the confinement is such that to inhibit the crystallization of water. Under 
this conditions, the dehydration of hydrophobic surfaces is characterized by the 
appearance of heterogeneities and cavitation. 

• The hydrophobic confinement affects the thermodynamics of water by lowering the 
T and increasing the P of the liquid-gas phase transition and of the TMD line. It 
also affects the LDL-HDL phase transition in a way that is possibly more complex. 

All these results are potentially relevant in problems such as the protein 
denaturation or the protein aggregation. Works are in progress to underpin and build 
up a theory of water at interfaces that could help us to acquire a better understanding 
of these subjects. 
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